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Abstract 


Symbolic computation techniques are used to derive some closed form expressions for 
an analytic continuation of the Euler-Zagier zeta function evaluated at the negative 
integers as recently proposed in [1], This approach allows to compute explicitly some 
contiguity identities, recurrences on the depth of the zeta values and generating 
functions. 


1 Introduction 


The multiple zeta functions, first introduced by Euler and generalized by D. 
Zagier [2], appear in diverse areas such as quantum field theory [5] and knot 
theory [7]. These are defined by 


Cr (^1 ) • • • ! ) 


E 

0<ki<---<kr 


kT...k 


( 1 ) 


where {rij} are complex values, and ([T]) converges when the constraints 


k 

Re (rir) > 1, and ^Re (rir+i-j) > k, 2 < k < r, (2) 

i=i 
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are satisfied (see 0). Their values at integer points n = (ni,..., Ur) satisfying 
(EJ are called multiple zeta values. An equivalent definition of these values is 


Cr (^ 1 ! ■ ■ ■ ) ^r) 
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fcl>0,...,fcr>0 


1 

fcr (ki + k2r...(ki + --- + krr' 


The sum of the exponents ni + • —h is called the weight of the zeta value, 

and the number r of these exponents is called its depth. 


Following the result by Zhao [8] that the multiple zeta function has an analytic 
continuation to the whole space C'’, several authors have recently proposed 
different analytic continuations based on a variety of approaches: Akiyama et 
ah [3] used the Euler-Maclaurin summation formula and Matsumoto [1] the 
Mellin-Barnes integral formula. 


B. Sadaoui [T] provided recently such analytic continuation based on Raabe’s 
identity, which links the multiple integral 




dx 


n = 


[l,+oo)^ (xi + Oi) ^ (xi + Oi + X2 + 02) ^ • (2:1 + Ol + • • • + Xr + Or-) 


to the multiple zeta function 


Z(n,z) 


E 


1 

{ki + Zi) ^ (fci -\- Zi-\- 1x2-\- ^ 2 ) ^... {ki Zi kj. -\- Zy.) 



B. Sadaoui uses a classical inversion argument to obtain an analytic contin¬ 
uation of the multiple zeta function defined at negative integer arguments 
— n = (—ni,..., —Ur) . The argument uses the following three steps: 

- the integral Fk (n) is computed for values of rii,..., that satisfy the con¬ 
vergence conditions ([2]), 

- the values n are replaced by — n in this result: it is then shown that (— n) 

is a polynomial in the variable a, 

- the variables a = (oi, ..., a^) are replaced by {Bi, ..., Br ), and each Bernoulli 
symbol Bk satisfies the two evaluation rules: 

evaluation rule 1: each power B^ of the Bernoulli symbol Bk should be 
evaluated as 

Bl ^ (3) 

the p—th Bernoulli number 

evaluation rule 2: for two different symbols Bk and Bi, k ^ I, the product 
B^Bf is evaluated as 

BlBi ^ B,B„ ( 4 ) 
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the product of the Bernoulli numbers Bp and Bq.\i k = I, the first rule applies 
to give the evaluation 

KK ^ Bp^q. 

Example 1 . An example of depth 2, appearing in |T], is now computed using 
the rules above. The integral (ui,?7,2) is explicitly computed and, replacing 
(ni,U 2 ) by (-ni,-n 2 ) gives 


d"ai,a2 (-ni, -^ 2 ) 


1 

77,2 + 1 


ri. 2+1 ni+n 2 + 2 —fc 2 ^2 

E E E 


/ll2=0 /i=0 ^2=0 
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Then substituting the variables Ui and a 2 with the Bernoulli symbols Bi and 
B 2 gives 


C2(-ni,-772) 


1 

772 + 1 


n2+l ni+n2+2—fc2 k2 

E E E 


k2=0 /i=0 l2=0 


f )12 + l\ /ni+712+2—fc2\ (k2\ 

\ ^2 / V _ h _ / \h) jgpg/2 

72-1 + 772 + 2 -^2 


Using the evaluation rules ([3]) and (jl]) for the Bernoulli symbols, the multiple 
zeta value of depth 2 at (—77i, — 772 ) is 


C 2 (- 7 ^ 1 , - 7 ^ 2 ) 


1 

772 + 1 


712 + 1711+712+2—^2 k2 

E E E 


/C2 = 0 /l=0 ^2=0 


^712 + 1^ ^7li+7l2+2-A:2^ ^^2^ 

7^1 + 772 + 2 - A:2 


Bi^ Bi^. 


The general case is given in [H eq. (4.10)] as the (2r — 


1) —fold sum0 




Cr(-77i,...,-77^) = (-l)’' 
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Bi. ...B, 


where k 2 , ■ ■ ■ ,kr > 0, Ij < kj for 2 < j < r and h < n + r + k and 


77 = ^77j, 

i=i 


k = Y.^r 

J=2 


( 6 ) 


A symbolic expression for (|5]) is proposed here. This is used as a convenient tool 
to derive some specihc zeta values at negative integers, contiguity identities for 
the multiple zeta functions, recursions on their depth and generating functions. 

1. This corrects a typo in [U eq. (4.10)] 
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2 Main result 


Introduce first the symbols Ci^ 2 ,...,k defined recursively in terms of the Bernoulli 
symbols Bi,... ,Br as 


C 


n 

1 


ton 

^1 fon _ 

) '-' 1,2 ~ 
n 


(Cl + B 2 T 
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and C 


n 


(Cl,2,...,A: 


+ Bk+l) 
n 


with the symbolic computation rule: 


C— symbols rule: All symbols Ci_ 2 ,...,A: are expanded using the above identities 
to express them only in terms of Bk- The evaluation rules ([3]) and (jl]) for the 
Bernoulli symbols are then applied. 

Example 2. For example, 


nn\nn2 _ nni 

is evaluated as 


{C. + B^r _ 1 ^ 1 ^ (n2\ BT 


n2 




^2 V k 


^^2 ^ V ^ ^ 


1 ^ fn2\ 

) ni + k 


Bn2—k- 


The next result is given in terms of this notation. 

Theorem 2.1 The multiple zeta values ^ at the negative integers (—ni,..., —nr) 
are given by 

Cr {-ni ,..., -nr) = n (-1)’"'= C”^^ 


( 7 ) 


k=l 


Proof 2.2 The inner sum in m. in its Bernoulli symbols version, 


E 

/l 


^n + r — k\ (k 2 


^kr 


h 


/g 




can be summed to 


(1 + (1 + B2)"^ ... (1 + Br)'^’-. 

The classical identit'j^ for Bernoulli symbols B + 1 = —B, with h defined in 
([6]) reduces this to 

( 8 ) 


\k2 


2 . this identity can be deduced from the generating function 

exp (zB) = 


exp (z) — 1 
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It follows that 


Cr (-n) = 


■ 1 / 
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Summing first over k 2 gives 


Y, n 


(' 
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i=2 


(EI.jni + r-i + 1-Etj+ifci) 


Cr (-n) = 


(- 1 / 


{Ur + 1) 
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... B^/ n 


k-i 


=j+i 


/Ci 


/C3 5 • • • 


/.3 (EEi Ki + r - i + 1 - E".i+, fe.)' 


The result now follows by summing, in order, over the remaining indices. 


Observe that the reduction ([8]) performed in the proof allows to restate a 
simpler version of Sadaoui’s formula ([S]) as the more tractable (r — 1) —fold 
sum 


Cr ( ^1) • • • ) 




-nr)-{ 1) Ym r\ 11 '' . . 

k 2 ,...,kr [n + r-k^ j=2 [Y.i=j Ui + r -j + 1 


... Bi^ 
- 
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Observe moreover that the derivation of o is unchanged if the symbols 
B„...,Br are replaced by a generalization of the Bernoulli symbol B, namely 
the polynomial Bernoulli symbol B + z dehned by 

{B + zf = Br,{z), 

the Bernoulli polynomial of degree n. The same proof as above yields the next 
statement. 


Theorem 2.3 The analytic continuation of the zeta function as given in 
can be written as 

r 

Cr {-ni, ..., -Ur, Zi, . . . ,Zr) = Y[ {^1, ■ ■ ■ , Zi) (10) 

i=l 


with 


cn ^ (^1 + Bi)"' ^ Bn (^l) ^n ^ (^1 (^ l) + ^2 + ^ 2 )' 


n 


n 


n 


and 


^l,2,...,A:+l (^1’ ■ ■ ■ ! 2 ^A:+i) — 


(Cl,2,...,A: {Zl, . . . , Zk) + Bk+1 + Zk+l) 
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3 A general recursion formula on the depth 


The methods above are now used to produce a general recursion formula on 
the depth of the zeta function. 

Theorem 3.1 The multiple zeta functions satisfy the recursion rule 

\Tlr Tlr-\-\ 


Qr (-n; z) = 




+ 1 


E 

fc =0 


k 


( 1) Cr—1 ( ^1) ■ ■ ■ ) ^r—l Brir+l—k (^r 


Introducing the new zeta symbol with the evaluation rri/e[l] 

Zr = Cr {-ni, . . . , -Ur-l, -Ur - k] z) , 

this recursion rule can he written symbolically as 

( n, z) (1) H Cl ( ^rj Zr—lj 

Ur + I 

Proof 3.2 Start from / fl^) and expand the last term 


( 11 ) 


( 12 ) 


Cf;+)iz,,...,Zr) = 


{f-ll-lr-l (^l5 • • • ) ^r-l) + (^r)) 


rir+l 


Ur + 1 


by using the binomial formula to produce 


Cr ( ^1) • ■ ■ ) ^li ••• 1 ^r') 






X qE'Ci (^1, ■ ■ ■, ^.-i) i^r) 


i=l 

Ur^l—k 


Then identify 


(n (z.,.... z,)) cztPi (z„.... .'.-i) 


as 


2+n, 1+fc ^ ^ ^ ^ -n^_ 2 , -Ur-l - k] z) 

to obtain the desired result. 

Using the symbol Z, this identity can be rewritten as 

Cr {-ni, -Ur, Zi,...,Zr) = - {B - 

77.j._|_i 


3. note that / 1 
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and the initial value 


Cl {-n; z) 

provides the stated reeursion. 


{z + 
n + 1 


4 Contiguity identities 

The multiple zeta function at negative integer values satisfies contiguity iden¬ 
tities in the variables. Two of them are presented here. 

Theorem 4.1 The zeta funetion satisfies the contiguity identity 

Cr ( ^1) • • • ) ^1) • • • ) Z^ — 1 , -|- 1) ( Til, • • • ) ^ri Z\i ■ ■ ■ , Z ,— i, Zj-^ -|- ( 1) (^Zf 2 ,— i) 

Example 3. In the case of the zeta function of depth 2, 

C2 (-ni, -ria, Zi, Z2 + 1 ) = C2 (-^i, -ris, Zi, Z2) + (z2 - 

and the second term is expanded as 

(-1)”'*' E ( 7 ) (-"1 - ■ 


Proof 4.2 Expand 


Cr (-ni, -Hr] Zi,..., Zr-l, Zr + 1) = ( 2 : 1 ) . . . C”" (Zl, . . . , Zr-2) 

an I ’ ’ 


Ur + 1 

flr + l 


X 


^ ^ (^1) • • • ) ^r-l) ^n^+l-fc (2^r + 1) 


and use the identity on Bernoulli polynomials 

Brir+l-k {Zr + 1) = -Srir+l-fc i^r) + {^r — k + 1) ^ 

to produce the result. 

The corresponding result for a shift in the hrst variable admits a similar proof. 
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Theorem 4.3 The depth-2 zeta funetion satisfies the contiguity identities 


C2{-ni,-n2,Zi + l,Z2) 


C2{-ni,-n2,Zi,Z2) + 


- 1 ) 


ni+n2 


712 + 1 


-Zi^Bn2+l {Zl + Z 2 ) ■ 


5 A Generating Function 


The generating function of the zeta values at negative integers is dehned by 

...ni Ur 

Fr{Wi,...,Wr)= \ ^ Cr (-77l, • • • , -Ur) ■ (13) 

A recurrence for is presented below. The initial condition is given in terms 
of the generating function for Bernoulli numbers 


B uj 

Fg (w) = y 

t„n< e"-l 


Theorem 5.1 The generating function of the zeta values at negative integers 
satisfies the recurrence 

Fr {wi, ...,Wr) = — [Fr-l (tCl, . . .,Wr-l) - Fb {-Wr) Fr-1 («^1, • • • ,tCr-2, ^r-l + ^r)] 

with the initial value 


Fi (wi) =- 

Wi 


-wiBi 


- 1 


1 - Fb 


Wi 


Moreover, the representation of the shift operator as exp o / {w) = 

f {w-\-a) and Fi{w,z) = —— — l] give the recursion symbolically 


as 


so that 


Fr {wi, . . . ,Wr) = Fi { Wr, 


d 


dWr-l 


O Fr-l {wi,.. .,Wr-l) , 


Fr {Wi, ...,Wr)=Fi{ Wr, 


d 


dWr-l 


oFi Wr-l, 


d \ ( d \ 

o- • -oFi ( W 2 , -tt— I 0+1 (tt’i) 


dWr-2 


dwi 


Proof 5.2 Start from 


Fr {wi. 


.Wr = 


E 

ni 5 ...nr 




■ ■ - w': 


(- 1 ) 


niH- hrir 




+1 


i=i 


-WjCi, 


np. ■ ■ - Url 



and expand 




n / 1 \n+l 


{-WrT (-1) 


n=0 


n\ n + 1 


(Cl ,_i + ■•■’’- 1 +^’') - l) , 

’ ’ Wr ^ ' 


to deduce that {wi,..., Wr) is 


(nc...■ 


>r-2 


Wr 




V=i 

= —Fr-1 (wi, . . . , Wr-l) - —Fb {-Wr) Fr-1 {Wl, . . . , Wr-2, Wr-1 + Wr) 
Wr Wr 


6 Shuffle Identity 


Multiple zeta values at positive integers satisfy shujfle identities, such as 

C2 (^1, n 2 ) + C2 (n 2 , ni) + Ci {ni + us) = Ci (^1) Ci (^^2) • 

The analytic continuation technique used in [1] does not preserve this identity 
at negative integers, while others do (for example, see [6]). The following 
theorem gives the correction terms. 


Theorem 6.1 The zeta values at negative integers as defined in m satisfy 
the identity 


C2 (-^1, -n2)+C2 (-^2, -ni)+Ci {-ni - n2)-Ci {-Ui) Ci {-n2) 


(_l)ni + l^l!^2! 
(ni + n2 + 2)! 

(14) 


Bni+n2+2- 


Remark 1. When ni + n 2 is odd, Bn,^+n 2+2 = 0 so that the shuffle identity 
(1T4|1 holds for (2 {—ni, —U 2 ) as expected, since the depth—2 zeta function is 
holomorphic at these points. 


Proof 6.2 Let 6 {wi,W 2 ) = F 2 {wi, W 2 )+F 2 {w 2 , wfi+Fi {wi + W 2 )-Fi {wi) Fi {W 2 ) ■ 
An elementary calculation gives 


6 {wi,W2) 


— + — I coth f i coth f 

Wl W2 2 \2/ 2 \2/ 

Wl + W 2 
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The expansions 


1 1 


Wi 2 

now produce 


i,-^-Y. 


+00 

Wi 


- {2k + 2)! 


B2k+2 (knd 


k=0 


Wi + W2 W2 


l>0 




W2 


+ 00 n 

6{w,,W2) = - E (-1)' f..2k+l+Y.,-l-l , 




k,i=0 (2fc + 2)! 

Identifying the coefficient of in this series expansion gives the result. 


7 Specific multiple zeta values 

This final section gives some examples of the evaluation at negative integers 
of the zeta function, obtained from ([5]) and ([12]). 

1: for depth r = 2, 


and 


C 2 ^ 

C 2 (0, -n) 


- 1 ) 

- 1 ) 


5, 


72+2 


1 B. 


n+l 


n + 1 


n + 2 2 n + ij ’ 

n+l 

[Bn+1 + -Sn+ 2 ] •; 


2: for depth r = 3, 


(15) 

(16) 


Cs (-n,0,0) = 


(-1)^ 


Bn+3 2 Bn+2 ^ 2 By^j^i 


2 Ln + 3 n + 2 3n + l. 


and 


Ca (0,-n,0) = 


(- 1 ) 


n+l 
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{n + 1) (n + 2) 


Tj Bn+3 

-Dn+2- — + 2-— 

n+l n+2 


(17) 


(18) 


3: as a hnal example, the recursion rule flTT|) is used to compute the value 
C3(0,0,-2) as 


C3(0,0,-2) 


^ (^ 3 ^° - 3B^Zl + 3BZl - 

i (S 3 C 2 (0,0) - 3 S 2 C 2 (0, -1) + 35+2 (0, -2) - C 2 (0, -3)) 


1 
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